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With present and future observations becoming of higher and higher quality, it is timely and neces- 
sary to investigate the most significant theoretical uncertainties in the predictions of inflation. We 
show that our ignorance of the entire history of the Universe, including the physics of reheating after 
inflation, translates to considerable errors in observationally relevant parameters. Using the infla- 
tionary flow formalism, we estimate that for a spectral index n and tensor/scalar ratio r in the region 
favored by current observational constraints, the theoretical errors are of order An/ |n — 1| ~ 0.1 — 1 
and Ar/r ~ 0.1 — 1. These errors represent the dominant theoretical uncertainties in the predictions 
of inflation, and are generically of the order of or larger than the projected uncertainties in future 
precision measurements of the Cosmic Microwave Background. We also show that the lowest-order 
classification of models into small field, large field, and hybrid breaks down when higher order 
corrections to the dynamics are included. Models can flow from one region to another. 

PACS numbers: 98.80.Cq 



I. INTRODUCTION 

Inflation pj] has become the dominant paradigm for 
understanding the initial conditions for structure for- 
mation and for Cosmic Microwave Background (CMB) 
anisotropics. In the inflationary picture, primordial den- 
sity and gravitational-wave fluctuations are created from 
quantum fluctuations, "redshifted" beyond the horizon 
during an early period of superluminal expansion of 
the universe, then "frozen" 0, 0, B- Perturbations at 
the surface of last scattering are observable as temper- 
ature anisotropics in the CMB, as first detected by the 
Cosmic Background Explorer satellite 0, The lat- 
est and most impressive confirmation of the inflationary 
paradigm has been recently provided by data from the 
Wilkinson Microwave Anisotropy Probe (WMAP) satel- 
lite, which marks the beginning of the precision era of 
CMB measurements in space 0|- The WMAP collabo- 
ration has produced a full-sky map of the angular vari- 
ations of the CMB to unprecedented accuracy. WMAP 
data support the inflationary mechanism as the mecha- 
nism for the generation of super-horizon curvature fluc- 
tuations. 

The CMB contains a wealth of information about the 
properties of the spectrum of primeval density perturba- 
tions and present data already allow to extract relevant 
informations about the parameters of single-field mod- 
els of infiation i.e. models whose inflation is driven 
by one scalar fleld, the inflaton. The following param- 
eters have been identified as important for accurately 
computing the expected anisotropy and for discriminat- 
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ing among different inflationary models: the power-law 
indices of the scalar and tensor perturbations n and 
riT respectively and the tensor-to-scalar amplitude ra- 
tio r — 16{Pt/ Pr)- Present data are consistent with a 
scale-invariant spectrum of scalar perturbation (n — I) 
and with an amount of tensor perturbations such that 
r < 0.5 jSj]. However, future CMB experiments will al- 
low an accurate determination of the properties the scalar 
spectrum. The satellite-borne experiment Planck Q , the 
proposed high-resolution version of CMBpol [13 and a 
polarized bolometer array on the South Pole Telescope 
will allow a determination of the spectral index n 
with a standard deviation of about 0.007, 0.003 and 0.01, 
respectively the same time a positive detection 

of the tensor modes through the i?-mode of CMB po- 
larization (once foregrounds due to gravitational lensing 
from local sources have been properly treated) requires 
r > 10"'^ ^M- While this limit is below the expected sen- 
sitivity, a tensor to scalar ratio of r 0.01 is well within 
the reach of presently feasible CMB observations. The 
proposed Big Bang Observer satellite has the potential 
to probe r - 0.001 E^. 

With present and future observations reaching a higher 
and higher quality, it becomes timely and necessary to 
investigate the most significant uncertainties on the the- 
oretical side as far as infiationary predictions are con- 
cerned. In this paper we study the impact of our ig- 
norance about the precise location on the infiationary 
potential corresponding to the observed perturbations. 
This is quantified by the number of e-foldings N before 
the end of inflation at which our present Hubble scale 
equalled the Hubble scale during inflation, the so-called 
epoch of horizon-crossing. Indeed, the determination of 
the number of e-foldings requires the knowledge of the 
entire history of the Universe. The expression for N can 
be written as 
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iV ~ 60 + ^In {-jit) + i In (Trh/IO^^ GeV) - ^ ln7 , 

(1) 

where Trh is the reheating temperature, 7 is the ratio 
of the entropy per comoving volume today to that after 
reheating and quantifies any post-inflation entropy pro- 
duction and we have assumed that there is no significant 
drop in energy density during the last stages of inflation. 
The main uncertainties in the determination of the num- 
ber of e-foldings are caused by our ignorance about the 
last two terms. The reheating temperature Trh after in- 
flation may vary from the Grand Unified Theory scale 
~ 10^^ GeV to 1 MeV, the scale at which nucleosynthe- 
sis takes place. In this range, the corresponding shift of 
N is about 14. Furthermore, long-lived massive particles 
of mass of the order of the weak scale are ubiquitous in 
string-inspired models (they are generically dubbed mod- 
uli) and may dominate the energy density of the Universe 
after reheating, leading to a prolonged matter-dominated 
epoch followed by a large amount of entropy release at 
the time of moduli decay [15| ,7^1. The corresponding 
shift in the number of foldings can be as large as 10. One 
can even envisage extreme situations where the reduction 
of the energy scale during inflation is so signiflcant that 
the shift in N is as large as 70 [iGj . 

Given the fact that the values of the inflationary ob- 
servables n, nr and r are evaluated at the value of N 
corresponding to the moment when the present Hubble 
scale crossed outside the horizon during inflation and that 
such a value is affected by a non-negligible uncertainty, 
we immediately conclude that the predictions of the infla- 
tionary observables are affected by unavoidable theoret- 
ical errors. How can we quantify them? Are they larger 
or smaller than the expected accuracy of forthcoming ex- 
periments? We propose to use the method of "flow" to 
gain some insight. The flow equations provide the deriva- 
tives of the inflationary observables with respect to the 
number of e-foldings as a function of the observable them- 
selves at any order in the so-called slow-roll parameters 
[TtL llsL H^] . For instance, to lowest order in slow roll 
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where ^ {m%J IQt:^) {H' H'" / H^) , mpi being the 
Planck scale, H the Hubble rate during inflation and 
primes indicate differentiation with respect to the in- 
flaton fleld. From this set of equations we can easily 
quantify - within any given single-field model of infla- 
tion - how our ignorance on the precise value of the 
number of e-foldings, quantified by a shift in the num- 
ber of e-foldings AiV, is reflected in the predicted value 
of the observable quantities. The expected uncertainties 
are model-dependent, but roughly speaking, we expect a 
theoretical error of the magnitude Ar ^ {dr / dN)^,/S.N 



and An ~ {dn/ dN)i,lS.N , where the derivatives are 
evaluated at a reference number of e-folding, e.g. at 
A^* ~ 60 corresponding to Trh ^ 10^^ GeV and 7 ~ 1 
and AA^ = (A^ — N^,). We conclude that - generi- 
cally - the error in the tensor-to-scalar ratio is of order 
Ar/r ^ 0.1 — 1, and the error in the spectral index is of 
order An/ |n — 1| ~ 0.1 — 1 (for AA ~ 14). These errors 
are of the order of or larger than the accuracy expected 
from future experiments. 

The paper is organized as follows. In Sec. II we discuss 
single-field inflation and the relevant observables in more 
detail. In Sec. Ill we discuss the infiationary model space, 
and in Sec. IV we describe the flow technique to quantify 
the theoretical errors in the infiationary predictions. In 
Sec. V we present our results and offer some analytical 
explanations Finally, in Sec. VI we present our conclu- 
sions. 



II. SINGLE-FIELD INFLATION AND THE 
INFLATIONARY OBSERVABLES 

In this section we briefly review scalar field models of 
infiationary cosmology, and explain how we relate model 
parameters to observable quantities. Inflation, in its most 
general sense, can be defined to be a period of acceler- 
ating cosmological expansion during which the universe 
evolves toward homogeneity and fiatness. This accelera- 
tion is typically a result of the universe being dominated 
by vacuum energy, with an equation of state p ~ —p. 
Within this broad framework, many specific models for 
infiation have been proposed. We limit ourselves here 
to models with "normal" gravity {i.e., general relativity) 
and a single order parameter for the vacuum, described 
by a slowly rolling scalar field (j>, the inflaton. 

A scalar field in a cosmological background evolves 
with an equation of motion 



4> + 3H^ + V' (0) = 0. 



(3) 



The evolution of the scale factor is given by the scalar 
field dominated FRW equation, 
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(4) 



We have assumed a fiat Fricdmann-Robertson- Walker 
metric. 



g^i, = diag(l,-a^ 



(5) 



where a^(i) is the scale factor of the universe. Inflation is 
defined to be a period of accelerated expansion, a > 0. A 
powerful way of describing the dynamics of a scalar field- 
dominated cosmology is to express the Hubble parameter 
as a function of the field cj), H = H {(/)), which is consistent 
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provided (j) is monotonic in time. The equations of motion 
become El IH El : 
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These are completely equivalent to the second-order 
equation of motion in Eq. (j^J. The second of the above 
equations is referred to as the Hamilton- Jacobi equation, 
and can be written in the useful form 
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where e is defined to be 
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The physical meaning of e can be seen by expressing Eq. 
(gl) as 
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so that the condition for inflation (a/a) > is given by 
e < 1. The scale factor is given by 
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Hdt 
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where the number of e-folds N is 
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whose value has been discussed in the Introduction. 

We will frequently work within the context of the slow 
roll approximation which is the assumption that the evo- 
lution of the field is dominated by drag from the cosmo- 
logical expansion, so that ~ and 
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The equation of state of the scalar field is dominated by 
the potential, so that p ~ — p, and the expansion rate is 
approximately 



H 



I 8n 
3r7ip[ 



(13) 



The slow roll approximation is consistent if both the slope 
and curvature of the potential are small, V, V" ^ V. 
In this case the parameter e can be expressed in terms of 
the potential as 
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We will also define a second "slow roll parameter" rj by: 
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Slow roll is then a consistent approximation for e, 77 <C 1. 

Inflation models not only explain the large-scale homo- 
geneity of the universe, but also provide a mechanism for 
explaining the observed level of inhomogeneity as well. 
During inflation, quantum fluctuations on small scales 
are quickly redshifted to scales much larger than the hori- 
zon size, where they are "frozen" as perturbations in the 
background metric. The metric perturbations created 
during inflation are of two types: scalar, or curvature per- 
turbations, which couple to the stress-energy of matter 
in the universe and form the "seeds" for structure for- 
mation, and tensor, or gravitational wave perturbations, 
which do not couple to matter. Both scalar and ten- 
sor perturbations contribute to CMB anisotropy. Scalar 
fluctuations can also be interpreted as fluctuations in the 
density of the matter in the universe. Scalar fluctuations 
can be quantitatively characterized by the comoving cur- 
vature perturbation P-ji . As long as the equation of state 
e is slowly varying, the curvature perturbation can be 
shown to be jlj 
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The fluctuation power spectrum is in general a function 
of wavenumber fc, and is evaluated when a given mode 
crosses outside the horizon during inflation, k = aH. 
Outside the horizon, modes do not evolve, so the ampli- 
tude of the mode when it crosses back inside the horizon 
during a later radiation- or matter-dominated epoch is 
just its value when it left the horizon during inflation. 
Instead of specifying the fluctuation amplitude directly 
as a function of k, it is convenient to specify it as a func- 
tion of the number of e-folds before the end of inflation 
at which a mode crossed outside the horizon. 
The spectral index n for Pji is defined by 
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din fc 

so that a scale-invariant spectrum, in which modes have 
constant amplitude at horizon crossing, is characterized 
by n = 1. 

The power spectrum of tensor fluctuation modes is 
given by 
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The ratio of tensor-to-scalar modes is then 
so that tensor modes are negligible for e <C 1. 
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III. THE INFLATIONARY MODEL SPACE 

To summarize the results of the previous section, infla- 
tion generates scalar (density) and tensor (gravity wave) 
fluctuations which are generally well approximated by 
power laws: 



P-R. (fc) cx k 



Pt (fc) cx fc"^ 



(20) 



In the limit of slow roll, the spectral indices n and hit 
vary slowly or not at all with scale. We can write the 
spectral indices n and to lowest order in terms of the 
slow roll parameters e and -q as: 
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(21) 

The tensor/scalar ratio is frequently expressed as a quan- 
tity r which is conventionally normalized as 
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The tensor spectral index is not an independent param- 
eter, but is proportional to the tensor/scalar ratio, given 
to lowest order in slow roll by 



-2e 



(23) 



This is known as the consistency relation for inflation. A 
given inflation model can therefore be described to lowest 
order in slow roll by three independent parameters, P-jz, 
Pt, and n. If we wish to include higher-order effects, we 
have a fourth parameter describing the running of the 
scalar spectral index, dn/dlnk. 

Calculating the CMB fluctuations from a particular 
inflationary model reduces to the following basic steps: 
(1) from the potential, calculate e and 77. (2) From e, 
calculate as a function of the field ip. (3) Invert N (0) 
to find (pN- (4) Calculate P-ji, n, and Pt as functions of 
(j), and evaluate them at (jj = cpM- For the remainder of 
the paper, all parameters are assumed to be evaluated at 

(j) = 4>N- 

Even restricting ourselves to a simple single-field infla- 
tion scenario^he number of models available to choose 
from is large 4] . It is convenient to define a general clas- 
sification scheme, or "zoology" for models of inflation. 
We divide models into three general types: large- field, 
small-field, and hybrid, with a fourth classification, lin- 
ear models, serving as a boundary between large- and 
small-field. A generic single-field potential can be char- 
acterized by two independent mass scales: a "height" A'', 
corresponding to the vacuum energy density during in- 
flation, and a "width" ji, corresponding to the change in 
the field value during inflation: 



y(0) 



(24) 



Different models have different forms for the function /. 
The height A is fixed by normalization, so the only free 
parameter is the width ^. 



With the normalization fixed, the relevant parameter 
space for distinguishing between inflation models to low- 
est order in slow roll is then the r — n plane. (To next or- 
der in slow-roll parameters, one must introduce the run- 
ning of n.) Different classes of models are distinguished 
by the value of the second derivative of the potential, 
or, equivalently, by the relationship between the values 
of the slow-roll parameters e and 77. Each class of mod- 
els has a different relationship between r and n. For a 
more detailed discussion of these relations, the reader is 
referred to Refs. ^4, 25]. 

First order in e and 77 is sufficiently accurate for the 
purposes of this Section, and for the remainder of this 
Section we will only work to first order. The generaliza- 
tion to higher order in slow roll will be discussed in the 
following. 



A. Large-Held models: —t<ri<t 

Large-field models have inflaton potentials typical of 
"chaotic" inflation scenarios 26] , in which the scalar field 
is displaced from the minimum of the potential by an 
amount usually of order the Planck mass. Such models 
are characterized by V" {(j)) > 0, and —e<ri<e. The 
generic large-field potentials we consider are polynomial 
potentials V {(j)) = {(f)/ fif , and exponential potentials, 
F(0)=A4exp(0//.). 

For the case of an exponential potential, V (0) cx 
exp(0//j,), the tensor/scalar ratio r is simply related to 
the spectral index as 
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(25) 



but the slow roll parameters are constant (there is no 
dependence upon N) and therefore no intrinsic errors of 
the observables n and r are expected in such a case. 

For infiation with a polynomial potential, V (</)) cx (/)^, 
we have 



71 - 1 = - 
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p + 2 



(26) 



so that tensor modes are large for significantly tilted spec- 
tra. By shifting the number of e-foldings by AiV one 
therefore expects 



A(77- 1) 

n - 1 



Ar 
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(27) 



From these relations we deduce that sizeable correlated 
theoretical errors should be expected for those large-field 
models characterized by large deviations from a flat spec- 
trum and by large values of the tensor-to-scalar ampli- 
tude ratio. Furthermore these errors increase with the 
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potential of the polynomial p. Of course, these state- 
ments are based on relations valid only at first order in 
the slow roll parameters. This means that for very large 
values of {n — 1) and r higher order corrections become 
relevant and may significantly alter the simple relations 



B. Small-field models: rj < ~e 

Small-field models are the type of potentials that arise 
naturally from spontaneous symmetry breaking (such as 
the original models of "new" inflation [27112^ ') and from 
pseudo Nambu-Goldstone modes (natural inflation [2^'). 
The field starts from near an unstable equilibrium (taken 
to be at the origin) and rolls down the potential to a 
stable minimum. Small-field models are characterized by 
V" {(j)) < and r] < —e. Typically e (and hence the tensor 
amplitude) is close to zero in small-field models. The 
generic small-field potentials we consider are of the form 
V {(p) = A'' [1 — ((/)/ /i)^] , which can be viewed as a lowest- 
order Taylor expansion of an arbitrary potential about 
the origin. The cases p = 2 and p > 2 have very different 
behavior. For p = 2, n— 1 ~ — (l/27r)(mpi//x)^ and there 
is no dependence upon the number of e-foldings. On the 
other hand 



C. Hybrid models: < e < ?7 



r = 8(f - n)exp[-l - 7V(f - n)] , (28) 



leading to 
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For p > 2, the scalar spectral index is 
n ~ 1 



N \p-2 



(29) 
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independent of (mpi/yi). Assuming fi < mpi results in an 
upper bound on r of 



r < 
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The corresponding theoretical errors read 
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(32) 



Due to the tiny predicted values of r, for small field mod- 
els one expects generically tiny errors in the tensor-to- 
scalar amplitude ratio, but sizeable errors in the spectral 
index. 



The hybrid scenario [30|, |31|, |33| frequently appears in 
models which incorporate inflation into supersymmetry. 
In a typical hybrid inflation model, the scalar field re- 
sponsible for inflation evolves toward a minimum with 
nonzero vacuum energy. The end of inflation arises as a 
result of instability in a second field. Such models are 
characterized by V" (</>) > and < e < ?/. We con- 
sider generic potentials for hybrid inflation of the form 
V {(j)) = A^[l + i(t)/^if] . The field value at the end of 
inflation is determined by some other physics, so there is 
a second free parameter characterizing the models. Be- 
cause of this extra freedom, hybrid models fill a broad 
region in the r — n plane. For {(jy^/fi) 3> 1 (where (f>N is 
the value of the inflaton field when there are e-foldings till 
the end of inflation) one recovers the same results of the 
large field models. On the contrary, when (0jv//i) ^ 1, 
the dynamics are analogous to small-field models, except 
that the field is evolving toward, rather than away from, 
a dynamical fixed point. This distinction is important to 
the discussion here because near the fixed point the pa- 
rameters r and n become independent of the number of 
e-folds iV, and the corresponding theoretical uncertain- 
ties due to the uncertainty in N vanish. However, there 
is an additional degree of freedom not present in other 
models due to the presence of the additional parameter 
(j)c. Therefore the theoretical uncertainties in the predic- 
tions of a generic hybrid inflation model are decoupled 
from the physics of reheating, and we do not consider 
such models further here. The distinguishing observa- 
tional feature of many hybrid models is 77 > and a blue 
scalar spectral index, n > 1. 

Notice that at first order in the slow roll parameters, 
there is no overlap in the r — n plane between hybrid 
inflation and other models. However, as we will explicitly 
show, this feature is lost going beyond first order: by 
changing N models can flow from the hybrid regions to 
other model regions; this feature is generic, models can 
flow from one region to another. Therefore it is important 
to distinguish between models labeled "hybrid" in the 
sense of evolution toward a late-time asymptote and the 
region labeled "hybrid" in the zoo plot. The lowest-order 
correspondence does not always survive to higher order 
in slow roll. 



D. Linear models: 77 = 

Linear models, V {(f) oc (/), live on the boundary be- 
tween large-field and small-field models, with V" {(j)) = 
and 77 — — e. The spectral index and tensor/scalar ratio 
are related as: 



r — — (1 — n) . 
3 ^ ' 



(33) 



For linear models, Eq. (|27|) applies. 
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This enumeration of models is certainly not exhaus- 
tive. There are a number of single-field models that do 
not fit well into this scheme, for example logarithmic po- 
tentials V {(j)) = Vo[l + {Cg^/Sn) In (^/^f)] typical of su- 
persymmetry Q , where C counts the degrees of freedom 
coupled to the inflaton field and 5 is a coupling constant. 
For this kind of potentials, one gets n — 1 ~ — (1/7V) and 
r ~ {2Cg'^ /tt^){1/N) corresponding to 

A(n-l) ^ Ar ^ _AN 
71-1 r N ' ^ ' 

Because of the loop-factor suppression, one typically ex- 
pects tiny theoretical errors in r, but sizeable uncertain- 
ties in n — 1. 

Another example is potentials with negative powers of 
the scalar field V ((p) = Vb [1 -I- a (nipi/cj))^] used in in- 
termediate inflation p33l | and dynamical supersymmetric 
inflation [33 . l35l |. Both of these cases require an auxiliary 
field to end inflation and are more properly categorized as 
hybrid models, but fall into the small-field region of the 
r — n plane. The power spectrum is blue being the spec- 
tral index given by n - 1 ~ 2(p+ l/p + 2)(2/(iVtot- A^)), 
where A^tot is the total number of e-foldings; the param- 
eter r turns out to be proportional to (n— 1)2(p+i)/(p+2) ^ 
Therefore, 

A(n- 1) _ p + 2 Ar _ AN 
n-1 ^ 2ip+l) r ~ Ntot-N- ^ ' 

Uncertainties in the spectral index can be sizeable if A'tot 
is close to N, but the theoretical errors in r are expected 
to be suppressed for small r. 

The three classes categorized by the relationship be- 
tween the slow- roll parameters as — e < 77 < e (large- 
field), rj < —e (small-field, linear), and < e < ry (hy- 
brid), cover the entire r — n plane and are in that sense 
complete (at least at first order in the slow roll param- 
eters) Figure ^ US shows the r — n plane divided into 
regions representing the large field, small-field and hybrid 
cases. Figure 12 |23l shows a "zoo plot" of the particu- 
lar potentials considered here plotted on the n — logr 
plane, along with projected errors from forthcoming ex- 
periments. For a given choice of potential of the form 

F(0)=AV(^^), (36) 

the parameter A is generally fixed by CMB normaliza- 
tion, leaving the mass scale and the number of e-folds 
N as free parameters. For some choices of potential, for 
example V cx exp (0//i) or ^ cx 1 — ((f>/ , the spectral 
index n varies as a function of /i. These models there- 
fore appear for fixed A^ as lines on the zoo plot. The 
inclusion of the uncertainty in N results in a broadening 
of the line. For other choices of potential, for example 
V (X 1— {4>/ fi)P with p > 2, the spectral index is indepen- 
dent of and each choice of p describes a point on the 
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n 

FIG. 1: Regions in the r — 71 plane corresponding to "large 
field", "small field", and "hybrid" models. 

zoo plot for fixed A^. The uncertainty in A^ turns each 
of these points into lines, which smear together into a 
continuous region in Fig. |21 Note that even if we include 
all of these uncertainties, the different classes of potential 
do not have significant overlap on the zoo plot, and it is 
therefore possible to distinguish one from another obser- 
vationally. Furthermore, for a given choice of potential, 
the uncertainties in r and n arising from the uncertainty 
in N are generally strongly correlated. This correlation 
will be apparent in the flow analysis presented below. Fi- 
nally, for particular choices of potential such as the expo- 
nential potential, inflation formally continues forever and 
the uncertainty due to A^ vanishes altogether, so there is 
no "smearing" of the line on the zoo plot. 



IV. FLOW EQUATIONS 

In this section we describe the flow equations which 
are a useful tool to quantify the theoretical errors in the 
inflationary observables due to our ignorance about of 
the number of e-foldings. 

We have defined the slow roll parameters e and rj in 
terms of the Hubble parameter H {(j)) as 

These parameters are simply related to observables r ~ 
16e, and n — 1 ~ 4e — 277 to first order in slow roll. Taking 
higher derivatives of H with respect to the field, we can 
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the same sign as H' {4>): 



n 



FIG. 2: Predictions of different potentials, including errors 
due to the uncertainty in A'^, plotted in the region of n — log r 
plane favored by current observation. The model curves, 
from largest r to smallest, are: V cx exp((;/>//i) (blue, top), 
V oc 0*" (magenta), V oc (/> (black), V oc 1 — (0//i)^ (red), 
1/ oc 1 — {(f)/ {P > 2) (green, lowest). The horizontal lines 
labeled with Qgw are the expected sensitivities of different 
proposed configurations of the Big Bang Observer satellite 
p^ . The hatched error bars are the observational uncertain- 
ties for: (a) a cosmic- variance limited temperature-only CMB 
measurement to £ = 1500 (outer, red), the Planck Surveyor 
satellite (middle, blue), and a hypothetical CMBPol-like ex- 
periment with the same angular resolution as Planck but three 
times better sensitivity (inner, solid black) [25|. The central 
value for the error bars shown is arbitrary. 



define an infinite hierarchy of slow roll parameters [3f 
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Here we have chosen the parameter a = 2ri — ~ n 
to make comparison with observation convenient. 

For our purposes, it is convenient to use N as the mea- 
sure of time during inflation. As above, we take and 0e 
to be the time and field value at end of inflation. There- 
fore, N is defined as the number of e-folds before the end 
of inflation, and increases as one goes backward in time 
{dt>0^dN < 0): 



_d_ 

dN 



dint 



mpi r- d 



(39) 



where we have chosen the sign convention that y/e has 



TOpL H' 



(40) 



2^ H ' 

Then e itself can be expressed in terms of H and N simply 
as, 

1 dH 



H dN 



e. 



(41) 



Similarly, the evolution of the higher order parameters 
during inflation is determined by a set of "flow" equations 

dN - ^(^ + 2^)' 



da 

dN 

dN 



-5ecr - 126^ + 2 (^Ah) , 
- 1 



-cr+ (£- 2)( 



rAH)+^+^AH.(42) 



The derivative of a slow roll parameter at a given order 
is higher order in slow roll. At the lowest order, this set 
of equations properly expressed in terms of observables 
reproduce equations 1^)1. 

A boundary condition can be specified at any point in 
the inflationary evolution by selecting a set of parame- 
ters e, a, ^Ah, • ■ ■ f^or a given value of N. This is sufficient 
to specify a "path" in the inflationary parameter space 
that specifies the evolution of the observables in terms 
of the number of e-foldings. Taken to infinite order, this 
set of equations completely specifies how a shift in the 
number of e-foldings is reflected in a shift of the slow roll 
parameters and, therefore, of the observables. Further- 
more, such a quantification is exact, with no assumption 
of slow roll necessary. In practice, we must truncate the 
expansion at finite order by assuming that the ^Ah are 
all zero above some fixed value of £. 

Once we obtain a solution to the flow equations 
[€{N),a{N), ^Ah(A^)], we can calculate the predicted val- 
ues of the tensor/scalar ratio r, the spectral index n, and 
the "running" of the spectral index dn/dlnk and how 
they change upon shifting the number of e-foldings by 
AA^. To lowest order, the relationship between the slow 
roll parameters and the observables is especially simple: 
r = 16e, n — 1 = a, and dn/dlnk = 0. To second order 
in slow roll, the observables are given by [s^ [s^ , 



r = 16e [1 - C (cr + 2e)] , (43) 
for the tensor/scalar ratio, and 

n-1 = cr-(5 - 3C) e^-i (3 - 5C) fie+i (3 - C) (^Ar) 

(44) 

for the spectral index. The constant C = 4(ln2+7) — 5 = 
0.0814514, where 7 ~ 0.577 is Euler's constant. Deriva- 
tives with respect to wavenumber k can be expressed in 
terms of derivatives with respect to N as ;38|1 
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The scale dependence of n is then given by the simple 
expression 



dn 
dink 



1 



dn 



1-e dN 



(46) 



which can be evaluated by using Eq. H44|l and the flow 
equations. 

It is straightforward to use the flow equations to obtain 
lowest-order estimates of the expected theoretical errors 
Ar and An in the predictions for r and n by adopting 
the simple approximation 



where 



dr 
dN 



16- 



Ar 
An 

de 
'dN 



AA^ 

dN 



16e (CT + 2e)) 

r [(n- 1) + (r/8)] 



(47) 



(48) 



and 



dn aa , , „ 9 „ /9 , \ 

dN-rn — 5— 12.^ + 2 f Ah) 



da 
dN 



We then have estimates for the uncertainties Ar and An 
in terms of the uncertainty in the number of e-folds AN: 



Ar/r 




An/(n-l)o_l 




FIG. 3: Lowest order estimates H5()|l of the theoretical er- 
rors Ar/r (top) and An (bottom) as functions of n and r, 
assuming AA*' — 14. 



r 

An 

'i-l 



[(n-l) + (r/8)] (AiV) 



16 



32n - 1 



(AiV), 



where we have taken ^Ah — 0. Figure |3| shows the error 
estimates from Eq. (|50f) as a function of r and n. These 
estimates indicate that the theoretical errors in n and r 
can be substantial, depending on where the model lives 
in the r-n plane. If we take the region roughly favored by 
current observation, r < 0.1 and |n — 1| < 0.1, we have 
upper bounds on the errors of order 



Ar 

— < 1 

r 

An 

< 1, 



1 



we expect substantial absolute error in the spectral index 
only in models which also predict a relatively large run- 

(50) ning dn/dlnk. These estimates indicate that the theo- 
retical errors in the inflationary observables can be signif- 
icant compared to the expected accuracy of future obser- 
vational constraints. Note that for a single-parameter set 
of models such as the large-field case (|27|) , there exists an 
iV-independent relation between r and n. Therefore the 
errors in the parameters are highly correlated. Such mod- 
els, as we have seen, are falsifiable by observation. Such a 
simple relation between the observables will not exist for 
models described by a larger number of parameters. In 
the next section, we present the results of a Monte Carlo 
analysis which extends these estimates to higher order in 
slow roll, effectively increasing the dimensionality of the 

(51) parameter space describing the potentials. 



where these bounds are saturated for r~0.1, |n— 1|~ 
0.1. Note in particular that, for n — 1 ~ 0.1, the fractional 
error in r is largely independent of the value of r. The 
fractional errors Ar/r and An/ |n — 1| are first order in 
slow roll. The absolute error An is second order in slow 
roll, and since 



An cx 



dn 
dink' 



(52) 



V. MONTE CARLO ESTIMATE OF 
THEORETICAL ERRORS 

In Sec. II VI we derived an analytical estimate of the 
theoretical errors in the observables n and r to lowest 
order in slow roll of 



Ar 
r 



< 1 
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While higher-order analogs of Eq. (|5U|I are in principle 
possible to derive using the flow equations, a comprehen- 
sive investigation of the effect of higher-order terms in 
slow roll is best accomplished using numerical techniques. 
In this section, we discuss the results of using a Monte 
Carlo evaluation of the flow equations to determine the 
errors Ar/r and An/ |n — 1| for a large ensemble of in- 
flationary models. 

Monte Carlo evaluation of the flow equations, intro- 
duced in Ref. has become a standard technique for 
investigating the inflationary model space. The princi- 
ple is straightforward: since the flow equations H42|) are 
first order differential equations, the selection of a point 
in the slow roll parameter space {e,77,^AH,...} serves to 
completely specify the evolution of a particular model in 
the space of slow roll parameters. For a model specified 
in this way, there is a straightforward procedure for de- 
termining its observable predictions, that is, the values 
of r, n— 1, and dn/dlnk a fixed number N e-folds before 
the end of infiation. The algorithm for a single model is 
as follows: 

• Select a point in the parameter space e, 77, 'Ah- 

• Evolve forward in time {dN < 0) until either (a) 
inflation ends, or (b) the evolution reaches a late- 
time fixed point. 

• If the evolution reaches a late-time fixed point, cal- 
culate the observables r, n—1, and dn/dlnk at this 
point. 

• If inflation ends, evaluate the flow equations back- 
ward N e-folds from the end of inflation. Calculate 
the observable parameters at this point. 

The end of inflation is given by the condition e — 1. In 
principle, it is possible to carry out this program exactly, 
with no assumptions made about the convergence of the 
hierarchy of slow roll parameters. In practice, the series 
of flow equations must be truncated at some finite 
order and evaluated numerically. The calculations pre- 
sented here are performed to eighth order in slow roll. ^ 
In effect, we are expanding our model space from the set 
of single-parameter potentials considered in Sec. IIIII to 
consider potentials with eight free parameters describing 
their shape [soj . 

We wish to determine how the uncertainty in the total 
number of e-folds N translates into a uncertainties in 
the observable parameters r and n. For models which 
reach a late-time attractor r — 0, n > 1 in the flow 




0.8 0.9 1 1.1 

n 

FIG. 4: Flow from iV = 46 (blue squares) to iV = 60 (red 
triangles) for an ensemble of fifty models generated via flow 
Monte Carlo, plotted in the n - r plane. The path traced 
by the flow indicates the level of theoretical uncertainty in- 
duced by the uncertainty AA'^. The diagonal lines indicate the 
boundares between small-field and large-field (green, solid) 
and large- field and hybrid (black, dashed). Models can "shift" 
class from A'' = 46 to A'^ = 60. The error bar at top right 
shows projected 2a measurement uncertainties in n and r for 
the Planck satellite. The central value for the error bar shown 
is arbitrary. 



space, the answer is trivial: since the observables are 
evaluated at a fixed point in the flow space, the shift in 
the observables with N by definition vanishes, and the 
theoretical uncertainty is in this sense negligible.^ We 
therefore concentrate on models dubbed "nontrivial" in 
the language of Ref. JlT}, that is models for which the 
dynamics carry the evolution through e = 1 and inflation 
naturally ends after a finite number of e-folds. For a given 
solution to the flow equations, it is simple to evaluate 
the effect of moving along th e "p ath" in flow space from 
iV = 46 to = 60. (Ref. contains an interesting 
analytic analysis of the dynamics of paths in the space of 
flow parameters.) 

Figures 01 \E\ and El illustrate the effect of shifting the 
number of e-folds in the space of the observables r, n, and 
dn/d\nk. We see that the effect of the theoretical error 
in A^ on the values of the observables is substantial, and 
at least qualitatively consistent with our rough estimate 
(|50|l and with our previous discussion for large and small 
field models. For large field models and for moderate 
value of (1 — rt) and r, for which the first order approxi- 
mations hold, we see that the errors increase respectively 
with {l~n) and r. Moving towards larger values of these 
parameters implies a substantial role played by higher or- 



^ The reader is referred to Ref. |17| for a more detailed discussion 
of the procedure used to stochastically evaluate the flow equa- 
tions. 



^ A more realistic hybrid-type model displaying this dynamics will 
likely be more complex, since the field may not yet have settled 
into the attractor solution at the appropriate point for calculat- 
ing cosmological observables. Such a situation is highly model- 
dependent, and we do not consider it further here. 
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FIG. 5: Flow from = 46 (blue squares) to = 60 (red 
triangles) for an ensemble of fifty models generated via flow 
Monte Carlo, plotted in the n - logr plane, showing the be- 
havior for small r. 
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FIG. 7: Theoretical uncertainties An and Ar/r calculated 
for an ensemble of 10,000 models with 0.9 < neo < 1.1. and 
reo > 0.001. The models cluster strongly in the region Ar/r ~ 
0.1 — 1 and An/ |n — 1| ~ 0.1 — 1, consistent with the lowest- 
order estimate (I50t . 
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FIG. 6: Flow from iV = 46 (blue squares) to iV = 60 (red 
triangles) for an ensemble of fifty models generated via flow 
Monte Carlo, plotted in the n - dn/dlnk plane. 
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FIG. 8: The absolute uncertainty in the spectral index An 
plotted versus the running dn/dlnk. As expected, the uncer- 
tainty in the spectral index becomes large for models which 
predict significant running. 



der corrections and errors can be very sizable. For small 
field models, we observe large displacements along the 
(n — l)-axis, but small ones along the r-axis. 

We can also see from Fig. 0]that the lowest-order clas- 
sification of models into small field, large field, and hy- 
brid breaks down when higher order corrections to the 
dynamics are included. Models which fall into the large- 
field region at iV = 46 can evolve into the small-field 
region at = 60, a behavior which was noted in Ref. 

Q. 

To obtain a more quantitative understanding of the 
theoretical error in the observables induced by the uncer- 
tainty AiV, we generate an ensemble of models using the 
fiow equations and calculate the observables at = 60, 
denoted rgo and ngo, and at N = 46, denoted and 
ri46. We retain only models which lie close to the region 
observationally favored by WMAP, 0.9 < ngo < 1-1, and 
we retain only models with a non-negligible tensor am- 
plitude, rgo > 0.001. For each model generated by the 



Monte Carlo, we then assign uncertainties in r and n as: 



and 



Ar 
r 



Ar 



?160 - "46 



(54) 



"60 



1 



(55) 



Figure 13 shows the above uncertainties calculated for an 
ensemble of 10,000 models. We see that the estimates 
Ar/r ~ An/ |n — 1| ~ 1 are robust even when calculated 
to higher order. Figure [S] shows the uncertainty An plot- 
ted against the running of the spectral index dn/dhik, 
showing the expected strong correlation between the er- 
ror An and the running. 
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VI. CONCLUSIONS 

We have considered the impUcations of the uncertainty 
in the reheat temperature on the observable predictions 
of inflation. The most convenient parameters for discrim- 
inating among inflation models are the tensor/scalar ra- 
tio r and the scalar spectral index n. These parameters 
are simply related to the inflationary slow roll parame- 
ters e and t], and this correspondence can be generalized 
to higher order in the slow roll expansion through the in- 
flationary flow equations. The uncertainty in the reheat 
temperature corresponds to uncertainty in the number of 
e-folds of expansion N during the inflationary epoch, 

TV ~ 60 + iln {-ut) + i In (Trh/IO^^ GeV) - ^ ln7 , 

(56) 

where Trh is the reheating temperature, 7 is the ratio 
of the entropy per comoving volume today to that after 
reheating. If we assume 7 is negligible, the large un- 
certainty in the reheat temperature corresponds to an 
uncertainty in the number of e-folds AiV ^ 14. This can 
be related to a theoretical uncertainty in the observable 
parameters r and n to lowest order by using the flow 
relations 



Ar 

r 

An 

'i-l 



[(n-l) + (r/8)] (AiV), 



16 



32n - 1 



(AiV). 



(57) 



For r and n in the region favored by current observa- 
tional constraints, these errors can in principle be large, 
Ar/r ~ 0.1 - 1, and An/ \n-l\ ~ 0.1 - 1. We also 
analyze the expected theoretical uncertainty in the in- 
flationary observables by Monte Carlo evaluation of the 
inflationary flow equations to eighth order in slow roll, 
and find results consistent with the lowest-order estimate 
above, but with considerable scatter in the models. We 
have numerically checked that these errors increase ap- 



proximately linearly with AA'^ as suggested by the lowest- 
order estimate. The absolute error in the spectral index 
can be compared with the expected 2cr uncertainty in the 
spectral index from the Planck satellite. An ^ 0.01, and 
is seen to be typically of the same order. 

We conclude that the theoretical uncertainties in the 
inflationary observables are generically of the order of or 
larger than the projected uncertainties in future precision 
measurements of the Cosmic Microwave Background, and 
represent a significant challenge for the program of us- 
ing observation to distinguish among the many differ- 
ent candidate models for inflation in the early universe. 
While the dependence of the inflationary observables on 
the number of e-folds N is certainly well known (and 
was, for example, taken into account by WHK in Refs. 
[TtI ll^), it has not been emphasized as the dominant 
source of theoretical error in the predictions of inflation. 
The error induced by the uncertainty in the reheat tem- 
perature and/or in the amount of entropy release after 
inflation (and thus in N) is typically much larger than 
errors in the quantities r and n due to using the slow roll 
approximation to calculate the primordial power spec- 
trum, a subject which has received considerable attention 
in the literature [H IH El El El . The expected 
uncertainties for any particular choice of potential are 
model-dependent: it will certainly still be possible to 
rule out models of inflation with future precision data. 
However, it will be necessary to move beyond the sim- 
ple lowest-order description of the inflationary parameter 
space which has so far been good enough. 
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